Isoscaling as a measure of Symmetry Energy in the Lattice Gas Model 
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The energetic properties of nuclear clusters inside a low-density, finite-temperature medium are 
studied with a Lattice Gas Model including isospin dependence and Coulomb forces. Important 
deviations are observed respect to the Fisher approximation of an ideal gas of non-interacting clus- 
ters, but the global energetics can still be approximately expressed in terms of a simple modified 
energy-density functional. The multi-fragmentation regime appears dominated by combinatorial 
effects in this model, but the isoscaling of the largest fragment in low energy collisions appears a 
promising observable for the experimental measurement of the symmetry energy. 
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In the different stages of the formation of a neutron 
star, from the explosion of the supernova to the cool- 
ing of the proto-neutron star, the baryonic matter is 
at finite temperature and densities different from the 
normal equilibrium density of nuclei. To understand 
these phenomena, a considerable experimental effort is 
presently devoted to the determination of the nuclear 
energy-density fuctional as a function of temperature 
and density [J, [2, H, [3l- Many investigations are con- 
centrated on the isovector properties of the system, the 
so-called symmetry energy, which behavior is still largely 
unknown. In the mean-field theory, the symmetry-energy 
term can be approximately reproduced as a siniple poly- 
nomial function of the density Csym{p) oc p^^] finite 
temperature leads just to a fractional occupation of 
single-particle levels, and as such does not modify the 
functional behavior of the interaction energy. 

The problem is that at low density the mean-field the- 
ory can be severely incorrect 0. In the heavy- ion col- 
lisions which are used to probe the density dependence 
of the symmetry energy, the system under study is typ- 
ically inhomogeneous; it is therefore not clear whether, 
even at thermal equilibrium, the associated energy func- 
tional only depends on the global density as obtained in 
a mean-field based picture. 

A simple formula has also been proposed^] to ex- 
tract directly the symmetry energy from measured clus- 
ter properties obtained in the fragmentation of two sys- 
tems of charge Zi,Z2, mass Ai,A2 at the same temper- 
ature T: 
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where a is the so-called isoscaling parameter that can be 
measured from isotopic yields [3]. 

However eq.([T]) has been derived in the framework of 
macroscopic statistical models f?] , where many-body cor- 
relations are supposed to be entirely exhausted by clus- 
terisation, and it appears to be strongly affected by con- 
servation laws and combinatorial effects [1, [lo| . Moreover, 
the Csym coefficient appearing in eq. ([1]) should correspond 



the symmetry energy only in the T — > limit. 

To progress on these issues, it is interesting to consider 
a microscopic model simple enough to be exactly solvable 
through Monte-Carlo simulations without any mean-field 
or independent-cluster approximation. In this paper, we 
study the temperature and density dependence of the 
symmetry energy coefficient and its connection to exper- 
imental observables in a Lattice Gas Model. 

Let us consider a system composed of N neutral and 
Z charged particles of mass m = 939MeV occupying 
a cubic lattice of ^ = 8000 cells with four degrees of 
freedom : one discrete variable at for isospin {ai = ±1 for 
protons (neutrons), ai — if the site is unoccupied), and 
three continuous variables pi for the momentum. The 
Hamiltonian of the system follows: 
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where < i,j > are nearest neighbor cells, e^iaj is the cou- 
pling between nearest neighbor (en ~ e-i-i = 0, ei_i = 
5.5MeV), Ic{= lAAMeV/ fm) is the Coulomb couphng 
between protons, and is the distance between sites 
i and j. The lattice spacing tq — 1.8 fm has been 
chosen such that a full lattice occupation corresponds 
to the saturation density of symmetric nuclear matter 
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This model has been already shown to be able to giv e 
a qualitative description of nuclear fragmentation 
Moreover the simplified version with only one type of 
uncharged particles is well-known to be isomorphous to 
the Ising model[l3|, which makes LGM a paradigm of 
first and second order phase transitions in finite systems. 

Calculations are made in the isobar canonical ensem- 
ble, which has been shown to be the correct canonical 
ensemble to describe unbound systems in the vacuum 
(l^ . The partition function reads: 
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to the symmetry free-energy llj , which is equivalent to where the sum runs over all the possible realizations (n) 
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of the system, and i?^*^"^ is the global extension of the 
system for each partition (n) defined as: 



(4) 



The partition sum ([3]) is numerically sampled for each 
given value of temperature and pressure with standard 
Metropohs techniques [l5l |. 

With only one type of particles, this model is well- 
known to exhibit a first order transition and a critical 
point {Tc,Pc), analogous to the liquid-gas transition. The 
phase diagram of a finite system can be obtained within 
the isobar canonical ensemble from the bimodality of the 
order parameter distribution At each pressure. 



the total energy distribution, as well as the distribution 
of the size of the heaviest cluster produced in each event, 
present two peaks of the same height at a temperature 
value, which is recognized as the transition tempera- 
ture [l8|. At this point the fluctuation are maximum, 
as shown in the left part of Fig[l] for a representative 
pressure. The ensemble of these transition points give 
the transition lines which are shown in the right part of 
FigllJ We can see that adding a short-range isovector 
coupling and a long-range repulsive interaction does not 
qualitatively modify the phenomenology of the liquid-gas 
transition (continuous line and full symbol) . In particular 
both fluctuations peak at the same temperature (dotted 
line), showing that the fragmentation transition has a 
flnitc latent heat also for charged systems. 

But the phase diagram is considerably enriched respect 
to liquid-gas. Two extra transitions appear at lower tem- 
perature which are speciflc to the nuclear phenomenol- 
ogy: inside the dashed curve, the system is splitted, with- 
out any energy jump, into two dominant fragments of 
similar size, which can be defined as hot fission. This re- 
sult is close to the findings of ref. . It is interesting to 
remark that bimodal distributions of the heaviest frag- 
ment have been recently observed experimentally [20]. In 
the following, we will only consider the systems above 
the residue-fission coexistence line. 

The presence of phase transitions implies that Lattice 
Gas systems are strongly inhomogeneous and cluster- 
ized. In this situation, it is not clear whether the en- 
ergetics of the system can be described by a macroscopic 
parametrization depending only on the average density, 
as in the mean-field approximation. 

To explore this issue, we try a liquid-drop inspired 
macroscopic parametrization for the interaction energy 
of the system: 

i?f„f(<5,p) = {a,{p) + c:y„,{p)6') A 

+ M + c%.MS^) A^'^ + a^Z^ (5) 

Here 5 is the isospin asymmetry 5 = {N — Z)/A, T is the 
temperature and p = A/ {4/3ttR^) is an estimation of 
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FIG. 1: Left side: fluctuation of the total interacting energy 
(upper part) and of the size of the largest cluster (lower part) 
as a function of the temperature at P = O.SPc, for a system 
with 75 neutrons and 54 protons. Right side: phase diagram. 
The full line is the liquid (L) -gas (G) coexistence line and 
the point is the critical point. The symbols with horizontal 
error-bars give fission (F) -residue (L) transition points. The 
dashed line is to guide the eye. 



the average density of the system, where the mean cubic 
radius — (i?"^*^"^) from eq.Q is calculated excluding 
the monomers {A = 1). 

This liquid-drop parametrization (LD) uses five macro- 
scopic density dependent parameters to be fitted to the 
model: is associated to the volume energy, cor- 
responds to the surface energy, Csj,„j and Csj,„j give the 
bulk and surface part of the symmetry energy, and etc 
corresponds to the Coulomb interaction. 

Eleven different systems of mass number A = 150 and 
isospin ratio ranging from S = — 1/3 to i5 = 1/3 are 
simulated at the pressure P/Pc = 0.3 and temperatures 
ranging from T/Tj = 0.88 to T/Tt = 1.1. This leads 
to an average volume variation between p/ pt = 0.7 and 
p/pt = 1.8. For each simulation, only one event out of 
2 X 10"* is kept to minimize auto-correlations, and av- 
erages are taken over 10^ events after a thermalization 
stage of typically 6 x 10'' events. The average energy cal- 
culated from the Metropolis simulation is confronted to 
the best fit obtained from eq. ^ on the upper part of Fig- 
ure [21 With a global /Ndof = 5, we can consider eq.([5]) 
as a reasonable approximation to the exact energetics of 
the systems. This is a non-trivial result, considering that 
the quantity p entering eq.(l5]) is never equal to the local 
density of the system, which is strongly fiuctuating. 

The density evolution of the macroscopic coefficients 
extracted from the best fit is plotted on the bottom part 
of Figure [21 The arrows give the values obtained when 
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the same fitting procedure is applied to the systems in 
their ground state 12, T^. We observe a decrease of all 
parameters with decreasing density, but the effect on the 
bulk terms is more important than the effect on the sur- 
face terms. Surprisingly, the contribution of the surface 
term to the symmetry energy appears to be negligible. 
This is consistent with the findings of ref. ^] . 

This result is encouraging for the experimental effort 
of extracting the nuclear matter Cgym out of nuclear col- 
lisions: in the framework of this model, the low val- 
ues extracted [l[ 0, 0] cannot be due to trivial surface 
effects (ii,^. 

The other interesting point is that the whole temper- 
ature dependence of the energetics is entirely embedded 
in the density dependence of the macroscopic parame- 
ters exactly like in the mean- field theory [j^ ]32i] . This is 
another non-trivial result, because different partitioning 
of the system could be associated to the same average 
spatial extension at different temperatures. What Fig. [2] 
demonstrates is that , at least in the framework of clas- 
sical physics and in the considered temperature-density 
range, the density functional approach can be a very good 
approach even in the presence of high-order correlations. 
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FIG. 2: Left part: Best fit obtained from the macro- 
scopic parametrization eq.([5]) versus the exact average energy 
at various average volumes and temperatures evaluated by 
Metropolis. Right part: Evolution of the macroscopic coeffi- 
cients of eq.(l5]) with the average density. The arrows give the 
zero temperature values of the macroscopic coefficients. 



This result suggests that it may be possible to extract 
the density dependence of the symmetry energy even if 
the system is clusterized. 

In order to explore this possibility, we now turn to ex- 
amine the cluster distribution by looking at the isoscal- 
ing observable. It is empirically well known that the ra- 
tio i?2i(A^, Z) of isotope yields Yi{N, Z) measured in two 
reactions labeled (1,2) at the same incident energy but 
different in isospin has an exponential behavior [8| ac- 



R2i{N,Z) = 



Y2{N, Z) 
Yi{,N,Z) 



cx exp {aN + (3Z) (6) 



where a and (3 are the isoscaling parameters. It is rea- 
sonable to imagine that an observable like a which mea- 
sures the fragment isotopic content should be sensitive 
to the symmetry energy, and in particular to its den- 
sity dependence [l|. Going a step forward, one may hope 
to extract Csym directly from the measured isoscaling 
through eq.([T|). 

It has already been observed that isoscaling is well re- 
spected in the LGM in the isochore ensemble [2^, [26j . 
The same is true also at constant pressure, as shown for 
a representative case in the left part of Fig. [31 The a 
parameter is almost constant with the fragment charge, 
and decreases with increasing temperature, which corre- 
sponds to decreasing density in the isobar ensemble. 

The symmetry coefficient extracted from eq.(IT]) using 
the value of a averaged from Z = 2 to Z = 7 is plotted 
on Fig. [3l We can see that the resulting parameter is 
almost constant and completely disagrees with the sym- 
metry energy of the model, as already observed in ref. [26[ . 
An alternative approximate formula was derived in the 
fragmentation regime in refs. [2l|, [i^l : 



a(Z) 
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{Z^l < A >l) 
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where Z is the charge of the considered fragment and 
< A>i\s the mean mass of this fragment as obtained in 
system z = 1, 2. In this expression Csym corresponds to 
the fragment symmetry energy, which may differ from the 
one of the source, for instance because of different surface 
effects[22, 23|. Different values of Csym{Z) (dashed hnes 
in Fig [3]) are observed for the different clusters charges 
applying eq.([7]). Indeed a is almost independent of Z, 
while the isotopic content of the fragments does depend 
on the fragment size: this is the well known fractiona- 
tion phenomenon[4, 23, 27|. Since the resulting Csym{Z) 
are constant for the different densities, eq.([7]) appears 
also inadequate to reproduce the symmetry energy of the 
model. As discussed in ref. [l3| , the weak sensitivity of the 
a parameter to Csym is due to the fact that light cluster 
yields are dominated by combinatorial probabilities, and 
do not reflect the thermodynamic properties of the sys- 
tem. 

The mass distribution of LGM is dominated by a large 
percolating cluster which is the order parameter of the 
fragmentation transition, and contains most information 
on the thermodynamics [28j. One may then expect that 
the isotopic distribution of the heaviest cluster produced 
in each event may be more sensitive to the symmetry 
energy of the fragmenting system. 

The result of applying eq.([7]) to the largest cluster is 
plotted on Figure[2]as a continuous grey line. We observe 
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a much better agreement, except close to the transition 
temperature. This may be due to the fact that very huge 
energy fluctuations are observed at the transition temper- 
ature in this canonical model (see Fig[T]), while the energy 
functional eq.fO depends only on average quantities. We 
expect that a better agreement will be obtained with a 
parametrization of the energy functional as a function of 
excitation energy and density, instead than temperature 
and density. This will further allow direct comparison 
with experimental data and will be the object of future 
investigations [l8l |. 
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FIG. 3: Left part: isotopic ratio as a function of the neutron 
number at T/Tt = 0.95 for the systems {N = 75, Z = 75) and 
(A'^ = 91, Z — 59). Right part: bold black line: symmetry en- 
ergy of LGM. The other lines give different estimation of Caym 
from isoscaling. Bold grey line: symmetry energy from the 
analysis of the biggest fragment. Dashed grey lines: eq.([7]), 
dashed black line: eq.([T}. 

To conclude, we have presented in this paper a study 
on the fragment properties at finite temperature and low 
density in the framework of a simple exactly solvable 
model . 

We have shown that even in thermodynamic configu- 
rations close to a phase transition, where the system is 
highly dishomogeneous and clusterized, the exact aver- 
age energy can be well described as a simple functional 
of the overall average density of the system. We have 
especially focussed our interest on the possible measure- 
ment of the density dependence of the symmetry energy, 
which is a topic of strong current interest in the nuclear 
physics community. We have shown that, in the frame- 
work of this model, the evolution of the symmetry energy 
term with the temperature and/or density can be traced 
with the help of the isotopic distribution of the largest 
cluster produced in each fragmentation event. 

Such measurements are presently undertaken by differ- 
ent experimental groups with the MARS recoil separator 
at Texas A&M[29] and the VAMOS spectrometer cou- 



pled with the INDRA 4 - tt array at GANIL [30] , and 
different experiments in this line are planned with future 
RiB's facilities jSlf. Such data should allow an important 
advance in the understanding of the functional behavior 
of the nuclear symmetry energy at finite temperature. 
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The same would not be true if we were to calculate the 
average free energy / = (e) — Ts instead than the av- 
erage energy (e), which is needed as an input of macro- 
scopic statistical models. This latter quantity would de- 
pend specifically on the temperature because of the en- 
tropy term. 



